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Coordinate  transformations  are  derived  from  global  Minkowski  coordinates  to  the  Fermi  co¬ 
ordinates  of  an  observer  moving  in  a  circle  in  Minkowski  space-time.  The  metric  for  the  Fermi 
coordinates  is  calculated  directly  from  the  tensor  transformation  rule.  The  behavior  of  ideal  clocks 
is  examined  from  the  observer’s  reference  frame  using  the  Fermi  coordinates.  A  complicated  relation 
exists  between  Fermi  coordinate  time  and  proper  time  on  stationary  clocks  (in  the  Fermi  frame) 
and  between  proper  time  on  satellite  clocks  that  orbit  the  observer.  An  orbital  Sagnac-like  effect 
exists  for  portable  clocks  that  orbit  the  Fermi  coordinate  origin.  The  coordinate  speed  of  light  is 
isotropic  but  varies  with  Fermi  coordinate  position  and  time.  The  magnitudes  of  these  kinematic 
effects  are  computed  for  parameters  that  are  relevant  to  the  Global  Positioning  System  (GPS)  and 
are  found  to  be  small;  however,  for  future  high-accuracy  time  transfer  systems,  these  effects  may  be 
of  significant  magnitude. 


I.  INTRODUCTION 

During  the  past  several  decades  there  have  been  significant  advances  in  the  accuracy  of  satellite-based  navigation 
and  time  transfer  systems  [1].  The  most  recent  example  of  such  a  system  is  the  Global  Positioning  System  (GPS), 
which  is  run  by  the  U.  S.  Department  of  Defense  [2-4].  The  GPS  is  based  on  a  network  of  24  atomic  clocks  in  space, 
from  which  a  user  can  receive  signals  and  compute  his  position  and  time.  The  order  of  magnitude  of  the  accuracy  of 
the  GPS  currently  is  10  m  and  20  ns  [5].  Future  navigation  and  time  transfer  systems  are  being  contemplated  that 
will  have  several  orders  of  magnitude  better  position  accuracy  and  subnanosecond  time  accuracy  [6-8]. 

Over  satellite  to  ground  distances,  precise  measurements  must  be  interpreted  within  the  framework  of  a  relativistic 
theory  [9].  The  design  of  high-accuracy  navigation  and  time  transfer  systems  depends  on  detailed  modeling  of 
the  measurement  process.  The  comparison  of  experiment  with  theory  in  general  relativity  is  more  subtle  than 
in  special  relativity  or  other  areas  of  physics  for  two  reasons.  First,  due  to  the  possibility  of  making  coordinate 
gauge  transformations,  it  is  often  not  clear  what  is  a  measurable  quantity  and  what  is  an  artifact  of  the  coordinate 
system  [10,11].  Second,  all  measurements  must  inprinciple  be  interpreted  in  the  comoving  frame  of  the  experimental 
apparatus  [12]  .  In  the  case  of  the  GPS,  even  theoretical  relativistic  treatments  do  not  interpret  the  measured  observable 
in  the  comoving  frame  of  the  apparatus,  (see,  for  example,  Schwarze  et  al.  [13]). 

The  principal  reason  for  investigating  relativistic  effects  in  detail  is  to  improve  the  current  accuracy  of  GPS  and 
to  create  future  time  transfer  and  navigation  systems  that  have  several  orders  of  magnitude  better  accuracy.  At  the 
present  time,  it  is  wellknown  that  small  anomalies  exist  in  position  and  time  computed  from  GPS  data.  The  origin  of 
these  anomalies  is  not  understood.  In  particular,  GPS  time  transfer  data  from  the  U.  S.  Naval  Observatory  indicates 
that  GPS  time  is  periodic  with  respect  to  the  Master  Clock,  which  is  the  most  accurate  source  of  official  time  for  the 
U.  S.  Department  of  Defense  [14-16].  Furthermore,  other  anomalies  have  been  found  in  Air  Force  monitor  station  data 
that  are  not  understood  at  present  [17].  The  work  below  is  a  first  step  in  revisiting  from  a  first  principles  approach 
some  of  the  physics  that  has  gone  into  the  design  of  the  GPS. 

A  theoretical  framework  known  as  the  tetrad  formalism  has  been  developed  for  relating  measurements  to  relativistic 
theory  [12,18-20].  However,  due  to  a  historical  lack  of  accurate  measurements  over  satellite  to  ground  distances,  in 
most  cases  it  has  not  been  necessary  to  apply  the  tetrad  formalism  in  the  area  of  satellite-based  systems.  This 
situation  is  rapidly  changing  due  to  improved  measurement  technology  [9,23-25]. 

In  the  tetrad  formalism,  measurements  are  interpreted  as  projections  of  tensors  on  the  tetrad  basis  vectors.  Con¬ 
sequently,  these  projections  are  invariant  quantities  (scalars)  under  transformation  of  the  space-time  coordinates. 
However,  these  projections  do  depend  on  the  world  line  of  the  observer  and  the  choice  of  local  Cartesian  axes  used 
by  the  observer  [20,12].  The  need  for  the  tetrad  formalism  to  relate  experiment  to  theory,  as  well  as  the  problem  of 
measurable  quantities  in  general  relativity,  is  well  discussed  by  Synge  [12],  Soffel  [21],  and  Brumberg  [22]. 
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The  tetrad  formalism  was  initially  investigated  for  the  case  of  inertial  observers,  which  move  on  geodesics 
[18-20,26-28].  However,  many  users  of  geodetic  satellite  systems  are  terrestrially  based  observers,  or  axe  based 
in  non-inertial  platforms,  such  as  missiles.  The  theory  for  the  case  of  non-inertial  observers  has  been  investigated 
by  Synge  [12],  who  considered  the  case  of  non-rotating  observers  moving  along  a  time-like  world  line,  and  by  others 
[29-33],  who  considered  accelerated,  rotating  observers.  For  arbitrary  observer  motions,  the  effects  are  indeed  very 
coinplicated,  and  the  general  theory  gives  limited  physical  insight.  For  this  reason,  in  this  work,  I  work  out  in  detail 
a.  simple  model  problem.  I  treat  the  particular  case  of  an  observer  moving  in  a  circle  in  Minkowski  space.  I  use  the 
tetrad  formalism  [12]  to  obtain  the  transformation  from  an  inertial  frame  in  Minkowski  space  to  the  accelerated  frame 
of  reference  of  an  observer  moving  in  a  circle.  The  observer  Fermi- Walker  transports  his  tetrad  basis  vectors.  The 
case  presented  here  corresponds  to  an  observer  kept  in  a  circular  trajectory  by  rockets:  i.e.,  he  is  not  on  a  geodesic. 
The  explicit  results  presented  below  form  a  particular  case  of  the  general  results  obtained  by  Synge  [12],  specialized  to 
zero  gravitational  field  and  circular  motion.  The  situation  treated  here  also  corresponds  to  a  reference  frame  that  has 
a  constant  orientation  with  respect  to  the  distant  stars  [34],  such  as  an  Earth-centered  inertial  reference  frame  used 
in  GPS  satellite  computations.  The  results  are  given  here  as  expansions  in  small  velocity  of  the  observer,  compared 
with  speed  of  light.  However,  the  coordinate  transformation  equations  are  not  expanded  in  a  power  series  in  time, 
and  hence  they  are  valid  for  arbitrarily  long  times.  ’ 

This  work  explores  the  apparent  behavior  of  ideal  clocks  as  observed  from  a  Fermi  coordinate  frame  of  reference  of 
an  observer  moving  in  a  circle  in  flat  space-time.  The  Fermi  coordinate  frame  is  the  closest  type  of  non-rotating  frame 
to  that  used  in  experiments  and  in  GPS.  In  particular,  a  Fermi  coordinate  frame  is  a  non-rotating  frame  of  reference, 
where  Coriolis  forces  are  absent,  and  light  travels  essentially  along  a  straight  line  [12].  The  observer’s  Fermi  frame  in 
this  work  is  analogous  to  the  Earth-centered  inertial  (ECI)  frame  used  in  GPS.  The  ECI  frame  moves  with  the  Earth 
but  maintains  its  orientation  with  respect  to  the  distant  stars. 

This  work  neglects  all  gravitational  field  effects  that  are  significant  for  the  actual  behavior  of  GPS  clocks.  As  such, 
this  work  treats  only  the  kinematic  effects  associated  with  the  circular  motion  of  the  Fermi  coordinate  frame  about 
a  central  point  (analogous  to  the  sun’s  position).  Since  gravitational  fields  are  neglected  in  this  work,  the  observer’s 
tetrad  is  Fermi- Walker  transported  along  a  time-like  world  line  (a  circle  in  3-d  space).  In  the  case  of  GPS,  the  origin 
of  the  ECI  frame  follows  a  geodesic  path  (the  Earth’s  orbit).  Therefore,  the  observer’s  Fermi  frame  in  this  work  is  2m 
accelerated,  non-rotating  frame,  while  the  ECI  frame  is  in  free  fall.  As  remarked  above,  the  effects  on  clocks  are  in 
general  very  complicated,  even  without  the  gravitational  field.  For  this  reason,  I  work  out  in  detail  the  specific  case 
u  “  circular  motion  in  Minkowski  space-time,  and  I  apply  the  results  to  parameters  that  are  relevant  to 


II.  FERMI- WALKER  TRANSPORT  DIFFERENTIAL  EQUATIONS 


Consider  a  time-like  world  line  C  of  an  observer  moving  in  a  circle  in  3-dimensional  (3-d)  space  given  by 


a:°  =  Iq  + 

4-  a  cos  (uu/c) 
=  Xq  +  osin  (uu/c) 


(1) 

(2) 

(3) 

(4) 


where  Xq,  Xq,  Xq,  and  are  constants,  u  is  a  parameter  along  the  world  line,  and  c  is  the  speed  of  light  in  an  inertial 
frame  in  vacuum.  The  observer  travels  in  a  circle  in  the  x-y  plane  (see  Figure  1).  Along  his  world  line,  the  observer 
travels  with  a  four  velocity 


and  has  an  acceleration 


(5) 
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where  is  the  affine  connection.  From  the  normalization  of  the  4-velocity,  u'ui  =  -1, 1  obtain  the  relation  between 
the  arc  length,  s  =  cr,  where  r  is  the  proper  time,  and  the  parameter  u, 


(7) 
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where  7  =  (1  —  i/^)  The  parameter  v  is  the  dimensionless  velocity,  given  by 

_  auj 
c 


(8) 


I  use  the  convention  that  Roman  indices  take  the  values  i  =  0, 1,2,3  and  Greek  indices  take  values  a  =  1,2,3  (see 
the  Appendix  for  the  conventions  used  in  this  work). 

In  Minkowski  space,  the  affine  connection  components  all  vanish:  =  0.  Hence,  the  explicit  expressions  for  the 

4~velocity  and  acceleration  are 


(*(s)  =  7(l,-^sin(^), 


(9) 


w*  =  — a7^  ^0,cos(^s),sin(^5),0^  (10) 

As  the  observer  moves  on  the  time-fike  world  line  C,  he  carries  with  him  an  ideal  clock  and  three  gyroscopes.  I 
use  the  notation  x*  for  Minkowski  space  coordinates.  At  some  initial  coordinate  time  x®,  the  observer  is  at  point  Pq 
at  proper  time  r  =  0,  corresponding  to  s  =  0.  On  his  world  line,  the  observer  carries  with  him  three  orthonormaJ 
basis  vectors,  a|2)j  where  a  =  1,2,3,  labels  the  vectors  and  i  =  0, 1,2,3  labels  the  components  of  these  vectors  in 
Minkowski  coordinates.  These  vectors  form  the  basis  for  his  measurements  [20]  (see  Figure  2).  The  orientation  of 
each  basis  vector  is  held  fixed  with  respect  to  each  of  the  gyroscopes’  axis  of  rotation  [35].  The  fourth  basis  vector  is 
taken  to  be  the  observer  4- velocity,  =  v^.  The  four  unit  vectors  a|^^j  form  our  tetrad  [12],  an  orthonormal  set  at 
Po,  which  satisfies 


9ij  ^(o)  ^Ib)  ~ 


(11) 


where  the  Sij  =  1  if  i  =  j,  and  zero  otherwise. 

At  a  later  time  s  =  cr,  the  observer  is  at  a  point  P,  The  observer’s  orthonormal  set  of  basis  vectors  are  related  to 
his  tetrad  basis  at  Pq  by  Fermi- Walker  transport.  Fermi- Walker  transport  preserves  the  lengths  and  relative  angles 
of  the  transported  vectors.  Furthermore,  Fermi- Walker  transport  of  basis  vectors  along  the  observer’s  world  line 
corresponds  to  keeping  the  vectors  as  parallel  as  possible  by  keeping  their  orientation  constant  relative  to  the  spin 
axes  of  the  gyroscopes  [35].  For  a  given  vector,  with  contravariant  components  /*  in  Minkowski  space,  its  components 
at  P  are  related  to  its  components  at  Pq  by  the  Fermi- Walker  transport  differential  equations  [12] 

sp 

(12) 

where 

(13) 

When  we  use  Eq.  (13)  to  transport  a  vector  p  that  is  orthogonal  to  the  4- velocity,  u*/i  =  0,  the  second  term  in  Eq.  (13) 
does  not  contribute.  We  refer  to  transport  of  such  space-like  basis  vectors  as  Fermi  transport,  and 
For  an  arbitrary  vector  p  perpendicular  to  the  4- velocity,  the  explicit  form  of  the  Fermi  transport  differential  equations 
for  the  world  line  in  Eqs.  (l)-(4)  is 

-pi/  cos  (7^)  p  -  7^1/  sin  (7^)  p  (14) 

pi/^  cos  (7^)  sin  (7^)  p  -h  pi/^  sin^  (7O  p  (15) 

-pi/^  cos^  (7^)  p  -  pv^  cos  (7O  sin  (7^)  p  (16) 

0  (17) 

where  the  components  p  are  functions  of  ^  and  the  dimensionless  proper  time  is  given  by 


dp_ 

di 

di 

dp_ 

di 

di 
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(18) 


The  differential  equations  in  Eqs.  (14)-(17)  are  nothing  more  than  the  equations  that  describe  the  Thomas  precession 
of  an  electron’s  spin  vector  as  it  moves  in  a  circular  orbit  around  the  nucleus  [36].  The  solution  of  Eq.  (14)-(17)  is 
given  by  [35] 


=  —72/^  cos  (7^^  4-  a)  +  /3  (19) 

f^  =  A  cos  (7^)  cos  (7^^  -h  a)  H-  i47  sin  (7^)  sin  (7^^  +  a)  (20) 

P  =  A  sin  (7^)  cos  (7^^  +  a)  —  ^47  cos  (7^)  sin  (7^^  +  a)  (21) 

f  =  8  (22) 

where  A,  a,  5  and  /3  are  integration  constants. 


III.  DETERMINATION  OF  THE  TETRAD 

At  the  point  Pq,  the  tetrad  basis  vectors  satisfy  Eq.  (11)., These  equations  constitute  12  relations  (since  Eq.  (11)  is 
symmetric  in  a  and  0)  for  16  components  .  However,  the  zeroth  member  of  the  tetrad  is  the  4- velocity,  which  is 
assumed  known: 


A|o)  =u‘  =  7(1,u1,v^,v3)  (23) 

and  is  given  by  Eq.  (9).  Equation  (11)  can  be  rewritten  as  two  equations: 

'^(a) '^(0)  “0)  a  =  1,2,3  (24) 

fl  =  1,2,3  (25) 

Substituting  the  explicit  form  of  the  Minkowski  metric  (see  the  Appendix)  into  Eq.  (24),  gives  a  condition  on  the 
zeroth  components  of  a  tetrad  vector  if  the  spatial  components  are  known  (Greek  indices  take  values  a,«  =1,2,3): 

(26) 

Then  use  of  Eq.  (26)  to  eliminate  the  time  components  from  Eq.  (25)  leads  to  an  orthogonality  relation  for  only 
the  spatial  components  of  the  tetrad  vectors 

~  ^ot(3  (27) 

The  general  idea  is  to  solve  Eq.  (27)  for  the  spatial  components  of  the  tetrad  at  s  =  0,  and  then  to  substitute  the 

spatial  components  into  Eq.  (26)  to  obtain  the  time  components  of  each  tetrad  vector.  Thus,  I  obtain  the  tetrad 

vectors  at  the  initial  time  corresponding  to  s  =  0.  For  each  tetrad  vector,  these  components  are  then  used  as  initial 
conditions  in  Eq.  (22  )to  determine  the  tetrad  components  at  point  P  at  finite  s. 

The  exact  solution  of  Eq.  (27)  can  be  obtained  by  noting  that  for  v"  =  0,  the  spatial  components  of  the  tetrad 

vectors  are  orthonormal,  and  the  solution  is  given  by 

^{a)  “  (28) 

I  assume  that  the  tetrad  unit  vectors  are  oriented  approximately  parallel  to  the  observer’s  local  Cartesian  x,  y  and  z 
axes.  For  small  «  1,  the  solution  of  Eq.  (27)  can  be  obtained  as  a  triple  power  series  in  u®,  by  iteration,  starting 
with  Eq.  (28)  as  the  first  approximation.  Having  solved  for  the  tetrad  components  as  a  power  series,  it  is  easy  to 
guess  the  exact  solution  to  be 


=  "rSaffV^ 

(29) 

4a) 

(30) 

4)  = 

=  7  (l,u^,u^,u^) 

(31) 
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where  the  components  of  dimensionless  velocity  are  given  in  Eq.  (9). 

Equations  (29)-(31)  give  the  components  of  the  observer’s  tetrad  basis  vectors  at  point  Pq.  At  point  P,  the  arc 
length  is  s  >  0,  and  Eqs.  (29)-(31)  give  the  initial  components  for  each  tetrad  vector.  These  initial  components  are 
used  in  the  general  solution  of  a  Fermi- Walker  transported  vector,  given  in  Eq.  (22).  The  tetrad  components  at  point 
P  at  finite  s  are  given  by 

.'^(0)  =  {t,  -7i'sin(70,  7i/cos(7^),  0}  (32) 

-^(i)  =  [-1^  sin(72  ^), cos(70  003(72 ^)  -|-  7  sin(70  sin(72  ^), 003(7^  e)  sin(70  -  7  003(7 e)  sin(72 0}  (33) 

“^(2)  =  {7 003(72  ^),  -7  cos(72 ^)  sin(7^)  -j-  cos(7 ^)  sin(72  ^), 7  cos(7^)  003(72 ^)  +  sin(7^)  sin(72 ^), 0}  (34) 

•^(3)  =  {0>0)0,1}  (35) 

Equations  (32)-(35)  are  valid  for  all  times  s  because  they  are  exact  solutions  of  Eqs.  (14)-(17)  and  satisfy  the 
orthogonality  relations  in  Eq.  (11). 


IV.  FERMI  COORDINATES 

Fermi  coordinates  are  defined  by  the  geometric  construction  shown  in  Figure  2.  Every  event  P'  in  space-time  has 
coordinates  x*  in  the  global  inertial  coordinate  system.  According  to  the  observer  moving  in  circular  motion  along 
a  time-like  world  line,  the  same  event  has  the  Fermi  coordinates  a  =  0,1, 2, 3.  The  first  Fermi  coordinate, 
A  V  >  _  a,  IS  just  the  proper  time  (in  units  of  length)  associated  with  the  event  P'.  The  proper  time  for  P'  is  defined 
as  the^ue  of  me  length  s  in  space-time  at  which  a  space-like  geodesic  from  point  P  (which  defines  s)  passes  through 

event  P  .  The  definition  of  s  requires  the  tangent  vector  of  this  space-like  geodesic,  be  orthogonal  to  the  observer 
4“Velocity  at  P, 


fi*Vi\p  =  0  (36) 

For  the  simple  case  of  Minkowski  space,  the  geodesic  is  simply  a  straight  line  connecting  the  points  P  P' 

M’  =  iV(P'-x*(75))  (37) 

where  x‘(u)  and  x*'  are  the  coordinates  of  P  and  P',  respectively,  in  the  global  inertial  frame,  and  AT  is  a  normalization 
constant  that  makes  fj,'  a  unit  vector.  The  orthogonality  condition  in  Eq.  (36)  is 

5iiM*(s)A^(,)(s)  =  0  (38) 

and  gives  m  implicit  equation  for  s  for  a  given  P'.  This  orthogonality  condition  gives  the  first  Fermi  coordinate  of 
the  point  P' 


a:(°)  =  s 


(39) 


and  Eq.  (38)  gives  s  as  an  implicit  equation: 

7s  =  x® '  -  Xq  -I- 1/  J(x^'  —  Xq)  sin(^^^-^)  -  (x2'  -  Xq)  cos(^^^^)j  (40) 

In  the  limit  of  small  speeds  of  the  observer,  i/  0,  Eq.  (40)  gives  s  =  x°'  -  xg.  For  u  «  1,  Eq.  (40)  can  be  solved 
for  s  by  iteration,  resulting  in  a  power  series  in  1/: 


=  s  =  Ax®  -f-  Ax^  sin(^^  cos(^^  ^) 

c  c 


i'  -  -Ax^u^ 

iU 


^  cos(^^)  +  2  Ax2  sin(^^)  Ax^  cos(^^)  -  Ax^  sin(^i^) 


,  Ax®  w , 
c 


-t- 0(i/)'‘  (41) 


where  I  use  the  notation 
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Ai’  =  x*'-Xo  (42) 

For  future  reference,  I  label  the  coordinates  of  P,  P',  and  Pq  in  the  global  inertial  frame  by  P'  = 

P  =  {x^,x^,x^,x^),  and  Pq  =  (a:8,a:J,a;g,xg). 

The  contravariant  spatial  Fermi  coordinates,  X^°‘\  a  =  1,2,3,  are  defined  to  be  [12] 

=VAi“^  =5ii<T(s)/x‘(s)r,(“^)  A^'^)(s)  (43) 

where  I  used  the  Minkowski  definition  of  space-time  distance,  a,  along  the  space-like  geodesic  between  P  and  P' 

=  Pij (x' '  ”  '  -  x^)  (44) 

where  gij  is  the  Minkowski  metric  and  =  r]ij  is  an  invariant  matrix,  which  is  also  equal  to  the  Minkowski 

metric  (given  in  the  Appendix).  The  parameter  s  in  Eq.  (43)  is  a  function  of  P'  coordinates  x*',  so  s  must  be 
eliminated  using  Eq.  (41),  so  that  are  explicit  functions  of  the  global  inertial  coordinates  x*'.  I  have  calculated 
this  transformation  correct  to  fourth  order  in  i/.  However,  since  the  expressions  are  complicated,  I  write  them  here 
correct  only  to  third  order  in  u: 

v(i)  A  1  If*!  A  1  .2Ax^u.  .  o  .  .2Ax^u;J  9 

=  Ax^  -  a  cos( - )  -(-  7  Ax^  -  Ax^  cos( - )  -  Ax^  sm( - )  1/^ 

c  '  A  "  c  ^  c  ^ 


-h— Ax®  -Ax^  -h  a  sin(— — ^ +  0{u) 


,Ax®a;^  1 


^(2)  A  2  ■  y^x''uj.  1.9  .9  ,2Ax®a;^  ,  1  .2Ax'^lj  ^ 

=  Ax^  -  a  sin( - )  -f  -  Ax^  -h  Ax^  cos( - )  -  Ax^  sin( - )  i/^ 

c  4  [  c  c 

+^Ax®  Ax^  —  a  cos(— ^ -h  0(1/)"^ 

The  inverse  transformation  can  be  obtained  from  Eqs.  (45)-(47)  by  iteration,  resulting  in 
x®'  -  xg  =  x(0)  +  [^(2)  cos(^^^)  -  sin(5^^)l  u+-  +  Oiuf 

x» '  -  xj  =  X(1)  +  a  cos(^^^)  +  1  A:(i)  -  Xd)  cos(^^)  -  sin(^^^)' 

x^'-xl  =  Xd)  +  a  sin(^i^)  +  1  X^)  +  X^)  cos(^^^)  -  X^d  sin(^ii^)'  +  0{u)^ 


x^'-xg  =xd) 


V.  METRIC  IN  FERMI  COORDINATES 

The  line  element  in  the  Fermi  coordinate  system  of  the  observer  is  written  as 

ds^  =  -G(ij)  dXd)  dXd)  (52) 

where  G(jj)  are  the  metric  tensor  components  when  the  Xd)  are  used  as  coordinates.  A  direct  calculation  of  the 
metric  tensor  components  from  the  tensor  transformation  rule 

dx^  dx^ 

G(ij)  =  Qkl  'qxU) 

using  the  transformation  Eqs.  (45)-(47)  and  the  Minkowski  metric  for  g^j,  gives 


G(oo)  —  -(1  +  C)^ 
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and  G(ai3)  =  where  C  is  given  by 


c  = 


Ax^  a; Ax®  Ax^  .  a; Ax® 

1 - cos( - ) - sin( - ) 

a  c  a  c 


+  0{v^) 


(55) 


Using  the  transformation  Eqs.  (45)-(47)  to  express  ^  as  a  function  of  the  Fermi  coordinates,  I  obtain: 

C  =  -^  cos(^X<°^)  +  sin(^x(o))]  1/2  +  0(i/^)  (56) 

The  result  in  Eq.  (54)  agrees  with  Synge’s  general  theory  [12]  of  Fermi  coordinates  specialized  to  flat  space.  Synge 
shows  that 

C  =  «;*  Af )  (57) 


To  third  order  in  v,  Eq.  (56)  agrees  with  the  general  theory  given  by  Synge  [12].  The  quantities  are 

the  covariant  Fermi  coordinates  and  =  lu*  A^^^  are  the  components  of  the  observer’s  4-acceleration  in  the  Fermi 
coordinate  system.  The  4-acceleration,  w*,  i  =  0, 1, 2, 3,  is  related  to  the  ordinary  acceleration,  a^,  0=1,2, 3,  by 


7^  (fx* 
c2  dfi 


(58) 


Using  this  relation,  C  can  be  written  in  terms  of  the  global  inertial  coordinates  of  P  and  P'  as 

C  =  ^<5,;5(a:“'-x“(s))a^ 


(59) 


For  the  circular  motion  treated  here,  from  Eq.  (56),  the  three  spatial  components  of  the  observer’s  acceleration  in 
Fermi  coordinates  are  given  by 


^cos(^X^®^),sin(^X^®)),0^  (60) 

where  I  have  dropped  fourth-order  terms  in  v.  To  third  order  in  i/,  the  observer’s  acceleration  is  just  the  classical 
centripetal  acceleration  for  an  observer  moving  in  a  circle. 

The  quantity  ^  is  a  potential  (analogous  to  a  gravitational  potential)  that  determines  the  rate  of  proper  time 
with  respect  to  coordinate  time  in  Fermi  coordinates.  Proper  time  is  kept  by  an  ideal  clock.  Coordinate  time  depends 
on  the  definition  of  the  reference  frame  and  coordinates  used  within  that  frame.  The  quantity  C,  given  in  Eq.  (59), 
depends  on  acceleration  a^,  and  the  difference  in  coordinates,  x“'  —  x‘^(s),  from  the  spatial  origin  of  the  Fermi 
coordinates.  Away  from  the  spatial  origin  of  Fermi  coordinates,  the  acceleration  produces  an  effective  gravitational 
potential  field,  (?  C,  which  leads  to  a  change  in  the  relation  of  coordinate  time  to  proper  time,  through  the  metric 
components 


VI.  APPARENT  BEHAVIOR  OF  CLOCKS  IN  FERMI  COORDINATES 


The  significance  of  proper  time  is  that  it  is  the  quantity  that  ideal  clocks  keep,  that  it  may  be  related  to  time  on  real 
clocks,  and  that  it  is  closely  related  to  measurable  quantities.  The  significance  of  coordinate  time  is  that  it  enters  into 
the  theory.  In  order  to  analyze  experiments  with  clocks,  we  must  relate  what  is  measured,  i.e.,  proper  time  intervals 
between  space-time  events,  to  coordinate  positions  and  times  of  these  events. 

Consider  an  ideal  clock  whose  world  line  is  given  in  Fermi  coordinates  by  =  A'(*)(X^®^).  At  coordinate  times 
and  X^q  \  the  clock  is  at  spatial  positions  and  x]q  \  respectively.  The  elapsed  proper  time  between  these 
two  events,  Ar  =  As/c,  is  given  by  the  integral  along  the  world  line  of  the  clock 


As  = 


dx(®)  dx(o) ; 


1/2 


(61) 
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A.  Stationary  Clock 


For  the  case  of  a  stationary  clock,  located  at  constant  Fermi  coordinate  position  the  relation  between  proper 
time  r  and  coordinate  time,  is  given  by  taking  =  0  in  Eq.  (52).  The  rate  of  proper  time  with  respect  to 

coordinate  time  is  then  given  by 

^  =  l  +  C(XW,X(“))  (62) 

Equation  (62)  gives  the  well-known  result  that  the  rate  of  proper  time  depends  on  the  location  of  the  clock  [42].  For 
the  simple  case  of  circular  motion  of  the  observer  with  C  given  by  Eq.  (56),  substitution  in  Eq.  (62)  and  integration 
from  to  Xg  ^  gives 

A,  =  4»)  -  X<»>  -  {x»)  [.in  (^4"’)  -  sin  (^4')]  '  H  (1^^"’)  '  }  <“> 

where  I  have  dropped  terms  of  0{v^).  Equation  (63)  shows  that  there  is  a  complicated  relation  between  elapsed 
proper  time  on  a  clock  located  at  Fermi  coordinates  and  Fermi  coordinate  time  X(°K  For  example, 

a  clock  on  the  X^^^  axis  (with  X^^^  =  0)  has  a  periodic  time  with  respect  to  Fermi  coordinate  time,  given  by 

As  =  X^^  —  (au;/c)X(^)  sin(a)X£^/c).  The  initial  conditions  in  Eq.  (63)  determine  complicated  phase  relations 
between  proper  time  and  coordinate  time. 


B.  Clock  in  Circular  Motion 


Now  consider  an  ideal  clock  that  is  a  distance  b  from  the  origin  of  Fermi  coordinates  and  the  clock  executes  circular 
motion  in  Fermi  coordinates  in  the  plane.  I  take  the  world  line  for  this  motion  to  be  given  by 


XW  =  6cos 

X(2)=6sin  (^XWj 
x(®)  =  0 


(64) 

(65) 

(66) 

The  clock  moves  in  the  X(^^-X^^  plane,  which  coincides  with  the  plane  defined  by  the  circular  motion  of  the  origin  of 
Fermi  coordinates  (see  Figure  2) .  The  situation  considered  here  is  a  simplification  of  the  problem  of  a  GPS  satellite 
(with  a  clock)  orbiting  the  Earth,  which  is  itself  orbiting  the  sun.  (In  the  case  of  GPS,  the  satellites  are  in  orbits  that 
are  inclined  at  55°  with  respect  to  the  Earth’s  equator.)  However,  as  mentioned  previously,  in  the  case  of  the  Earth 
and  GPS  satellite,  both  travel  on  geodesics.  Since  we  are  neglecting  the  effects  of  gravity  in  the  problem  considered 
here,  the  Fermi  observer  and  orbiting  clock  are  both  Fermi- Walker  transported  (e.g.,  with  the  aid  of  rockets). 

Substitution  of  the  satellite  clock’s  world  line  into  Eq.  (61)  leads  to  the  following  relation  between  elapsed  proper 
time  observed  on  the  clock  and  Fermi  coordinate  time 


As  = 


1  /bn\^  (0)  (0).  ab  r.  (0)\  . 

-  2  (tJ  W’  -^V)  - [sm(— (— X 


The  rate  of  proper  time  with  respect  to  coordinate  time  is 

^  /a;\ 2 

—  a  6  cos 


ds 

d^ 


-li 


—  u 


X(o) 


(67) 


(68) 


The  proper  time  has  a  constant  rate  offset  from  coordinate  time,  represented  by  the  first  term  in  Eq.  (67).  In 
addition,  the  proper  time  has  a  periodic  component  with  respect  to  coordinate  time,  given  by  the  second  term  in 
Eq.  (67). 

To  demonstrate  the  magnitude  of  these  effects,  I  use  parameters  that  correspond  to  Earth-Sun  distance  and  GPS 
satellites.  I  compute  the  magnitudes  of  the  following  terms  with  the  values  shown  in  Table  I: 


1 

2 


=  8.348  X  10"“ 


(69) 
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(70) 


1 

2 


27r 

a 


=  3.597  X  10-® 


sec 


Equation  (69)  gives  the  constant  rate  offset  of  the  moving  clock  with  respect  to  Fermi  coordinate  time,  due  to  time 
dilation.  In  GPS,  this  term  leads  to  a  net  slowing  down  of  GPS  clocks  by  approximately  7  /xs  per  day  [37],  which 
results  from  multiplying  the  value  given  in  Eq.  (70)  for  a  single  revolution  by  a  factor  of  2  (since  the  GPS  satellites 
make  approximately  two  revolutions  per  Earth  day.) 

The  second  term  in  Eq.  (67)  is  a  periodicity  of  the  proper  time  with  respect  to  Fermi  coordinate  time  with  amplitude 


ab 

c  “  UJ 


=  12.05  X  10"®  s 


which  corresponds  to  an  amplitude  in  the  rate 


(71) 


=  1.755  xl0“^2g  (72) 

The  proper  time  given  by  Eq.  (67)  is  periodic  with  respect  to  Fermi  coordinate  time  with  the  difference  frequency, 
Q  —  uj,  which  is  the  difference  frequency  between  satellite  and  observer  rotations.  This  is  a  kinematic  effect  due  to 
the  type  of  coordinate  system  used  (Fermi  coordinates)  and  the  fact  that  the  observer  (Fermi  coordinate  frame)  is 
moving  along  an  arc  of  a  circle,  and  therefore  experiences  an  acceleration.  Note  that  the  periodic  term  vanishes  when 
the  observer  angular  velocity  a;  =  0.  If  the  observer  were  moving  in  a  straight  line  and  hence  had  zero  acceleration, 
=  0,  the  periodic  effect  would  also  be  absent,  as  can  be  seen  from  Eq.  (57). 


C.  Sagnac-Like  Effect  for  Portable  Clocks 


Sagnac  [38-40]  demonstrated  that  there  exists  a  phase  shift  between  two  counter  propagating  light  beams  on  a 
rotating  platform  and  that  this  phase  shift  depends  on  the  angular  frequency  of  rotation  of  the  platform.  In  the 
Fermi  coordinates  of  the  observer  moving  in  a  circle,  there  exists  a  Sagnac-like  effect  for  two  portable  clocks  that  are 
orbiting  in  opposite  directions. 

Consider  a  clock  that  moves  eastward  along  a  circle  of  radius  b  in  the  equatorial  plane  of  the  Fermi  coordinates 
with  world  line  given  by  Eqs.  (64)~(66)  (see  Figure  2).  The  speed  of  the  clock  in  Fermi  coordinates  is  bCt/c.  Assuming 

that  the  clock  begins  its  journey  at  point  A  at  time  =  0  and  ends  it  at  =  27rc/n,  which  corresponds  to  one 
revolution,  the  proper  time  that  elapses  on  the  clock  is  given  by  Eq.  (67)  as 


As^  = 


ab  uP" 
c  Cl  — u 


(73) 


The  first  term  in  Eq.  (73)  is  the  time  dilation  due  to  the  speed  bCl/c  of  the  clock.  The  second  term  depends  on  the 
path  traversed. 

Next,  consider  an  identical  clock  that  starts  at  point  A  at  time  =  0  but  moves  westward.  The  world  line  of 
this  clock  is  then  given  by  Eq.  (64)-(66)  but  with  Cl  =  -|n|.  The  proper  time  on  this  clock  is.given  by 


As^  = 


1 


ab  up 

—  TKi - sin 

c  n  -ho; 


(74) 


When  the  east-moving  and  west-moving  clocks  have  each  made  one  revolution  and  they  have  returned  to  their  starting 
positions  at  point  A,  the  difference  in  their  times  is  given  by  Ar^  -  Ar_  =  (As+  -  As^)/c.  For  the  case  applicable 
to  GPS,  b  «a  and  u;  <<  |f]|,  the  difference  of  proper  time  on  the  two  clocks  is  given  by 


=  2^ 
a  \  c 


V 


.  b  /aWN2  LJ 


(75) 

(76) 


where  in  the  last  line  I  have  dropped  terms  of  0{u}^/Cl^).  This  Sagnac-type  effect  is  due  to  the  non-inertial  nature  of 
the  Fermi  coordinate  system.  The  origin  of  this  effect  can  be  understood  from  the  point  of  view  of  the  global  inertial 
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coordinate  system.  If  the  Fermi  coordinate  origin  moved  along  a  straight  line  in  Minkowski  space  (rather  than  a  circle), 
then  there  would  be  no  acceleration  and  =  0.  In  this  case,  the  east-moving  and  west-moving  clocks  would  have 
zero  time  diflFerence  after  one  revolution.  However,  since  the  origin  of  the  Fermi  coordinates  moves  along  a  circle,  the 
acceleration  picks  out  a  direction  in  the  Fermi  coordinate  space-time,  and  the  symmetry  between  the  two  clock  world 
lines  is  broken.  Consequently,  the  elapsed  proper  time  is  different  on  the  two  clocks  when  they  are  compared  after 
one  revolution.  Note  that  the  time  difference  on  the  two  clocks  increases  with  increasing  angular  velocity  of  rotation 
cj,  analogous  to  the  optical  Sagnac  effect  [38^0].  This  increase  is  due  to  the  fact  that  the  breaking  in  symmetry  is 
due  to  the  magnitude  of  the  acceleration,  see  Eq.  (57).  Using  the  values  of  the  parameters  appropriate  to  GPS 
satellites,  given  in  Table  1, 1  find  the  magnitude  of  the  time  difference  on  the  two  clocks  after  one  revolution  to  be 

T+  —  r_  =  2.066  X  10”^^  s  (77) 

This  is  a  small  effect  that  vanishes  when  the  origin  of  Fermi  coordinates  moves  along  a  straight  line  or  when  cj  =  0. 


VII.  COORDINATE  SPEED  OF  LIGHT 

The  three-dimensional  space  of  Minkowski  space-time  is  homogeneous  and  isotropic,  and  consequently,  the  speed  of 
light  is  a  constant.  However,  for  the  case  of  Fermi  coordinates  of  an  observer  moving  in  a  circle,  the  symmetry  of  the 
3-d  space  is  broken  by  an  acceleration.  Consequently,  the  space  is  not  homogeneous  and  isotropic,  and  the  coordinate 
speed  of  light  depends  on  position  in  the  3-d  space. 

For  a  general  space-time  metric  gij ,  the  coordinate  speed  of  light  Vc  in  the  direction  of  a  unit  vector  is  given  by 


where 


Vc 


1  -b  TaC"  ^ 


lot  = 


(78) 


(79) 


and  the  unit  vector  satisfies  the  relation  7^/?  =  1,  where  the  3-d  spatial  metric  is  given  by  [42] 


loL0  —  “ 


gOg  90(3 
500 


(80) 


Equation  (78)  shows  that  for  a  general  space-time  metric  the  speed  of  light  depends  on  direction  e“  and  on  position 
Md  time,  through  the  metric  components.  A  non-zero  off-diagonal  metric  component  leads  to  an  anisotropy 
(directional  dependence)  of  the  coordinate  speed  of  light.  The  metric  for  Fermi  coordinates  given  in  Eq.  (54)  has 
vanishing  off-diagonal  terms,  7^  =  0,  and  consequently  the  coordinate  speed  of  light  is  isotropic.  However,  the  metric 
component  G(oo)  has  a  non-trivial  position  and  time  dependence,  and  Eq.  (78)  gives  the  coordinate  speed  of  light: 


Vc  =  c(H-C)  (81) 

where  C  is  a  function  of  Fermi  coordinate  position  and  time,  given  by  Eq.  (56).  Therefore,  I  can  write  the  coordinate 
dependence  of  the  speed  of  light  as 


T  =  -T  +  0(,.*)  (82) 

This  variable  speed  of  light  is  a  kinematic  effect  due  to  the  non-inertial  nature  of  the  Fermi  coordinate  system  Note 
that  at  any  given  position  X  =  the  speed  of  light  depends  on  time  in  a  periodic  way,  with  the 

orbital  period  of  the  reference  frame,  u.  This  expression  can  be  rewritten  as 

Ac  . 

-  =  «(-)  X.n-fO(i.^)  (83) 

where  the  unit  vector  n  =  (-cos(^X(°))  ,-sin(fX(o))  ,0)  is  in  the  observer’s  orbital  plane  and  points  in  the 
direction  of  acceleration,  toward  the  center  of  the  circle  (see  Figure  3).  As  the  Fermi  reference  frame  moves  around 
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the  circle,  the  vector  n  always  points  toward  the  circle’s  center  in  Fermi  coordinates.  At  any  given  time  Eq. 
(83)  shows  that  to  second  order  in  u,  in  the  plane  given  by 

X  •  n  =  0  (84) 

the  speed  of  light  does  not  differ  from  its  vacuum  inertial  frame  value  c.  In  the  3-d  half-space  containing  positions  X 
closer  to  the  center  of  the  circle,  the  speed  of  light  is  increased.  For  the  3-d  half-space  containing  positions  X  that 
are  further  from  the  center  of  the  circle,  the  speed  of  light  is  decreased.  As  time  increases,  the  dividing  plane, 
given  by  Eq.  (84),  rotates  at  angular  frequency  uj  around  the  circle,  with  the  origin  of  Fermi  coordinates  remaining 
the  point  at  which  the  plane  is  tangent  to  the  circle. 

The  magnitude  of  the  inhomogeneous  variation  in  the  speed  of  light  in  Eq.  (82)  is  small.  For  the  kinematics  of 
GPS,  when  I  use  the  numbers  given  in  Table  1,  with  the  coordinates  X^^^  ^  b,  so  the  effect  has  a  magnitude 

w2 

—  6  =  1.755  x  (85) 

a 

and  has  a  complicated  position  and  time  dependence,  given  by  Eq.  (82).  Recently,  GPS  data  have  been  used  to 
obtain  an  upper  bound  on  the  anisotropy  of  the  speed  of  light:  Ac/c  <  10"®  (see  Wolf  and  Petit  [43]  and  references 
contained  therein).  Given  the  present  accuracy  of  GPS,  the  small  inhomogeneity  of  the  speed  of  light  given  by  Eq. 
(82),  is  probably  not  measurable. 


VIII.  SUMMARY 

I  have  considered  the  problem  of  an  observer  moving  in  a  circle  around  a  central  point  in  flat  space  (no  gravitational 
fields).  The  observer  carries  with  him  three  gyroscopes  and  uses  them  to  define  a  set  of  non-rotating  axes.  He  keeps 
his  axes  at  a  constant  orientation  with  respect  to  the  gyroscopes’  spin  axes.  The  observer’s  axes  are  Fermi  transported 
along  his  world  line  and  Fermi  coordinates  are  defined.  The  Fermi  time  coordinate  is  just  the  observer’s  proper  time, 
according  to  the  construction  shown  in  Figure  1.  Fermi  coordinates  play  the  role  of  ECI  frame  coordinates  in  GPS.  I 
investigate  the  kinematic  effects  seen  on  ideal  clocks  from  the  point  of  view  of  these  Fermi  coordinates.  Specifically, 
I  consider  two  cases:  stationary  clocks  (with  respect  to  Fermi  coordinates)  and  satellite  clocks  in  circular  orbit  with 
respect  to  Fermi  coordinates.  The  proper  time  interval  between  two  events  on  the  stationary  clock,  given  by  Eq. 
(63),  is  a  complicated  periodic  function  of  the  clock’s  Fermi  coordinate  time  and  coordinate  position  in  the  plane 
perpendicular  to  the  rotation  axis.  I  considered  next  a  clock  that  orbits  the  Fermi  coordinate  origin  (similar  to  a 
satellite)  in  the  equatorial  plane  of  rotation  and  computed  the  relation  of  proper  time  to  coordinate  time.  I  find  that 
the  proper  time  on  the  clock  has  time  dilation  (constant  rate  offset)  plus  a  periodic  dependence  on  Fermi  coordinate 
time  at  the  difference  frequency,  —  cj,  between  satellite  clock  orbital  frequency,  f2,  and  Fermi  coordinate  frame 
orbital  frequency,  a;. 

Next,  I  looked  at  the  time  difference  of  two  counter-orbiting  clocks  in  the  equatorial  plane  of  rotation,  one  west- 
moving  and  the  other  east-moving.  After  one  revolution,  the  difference  in  proper  time  on  the  clocks  shows  a  small 
Sagnac-like  effect,  increasing  with  the  angular  frequency  of  revolution  of  the  Fermi  frame,  a;,  given  by  Eq.  (76)  . 

Finally,  I  computed  the  coordinate  speed  of  light  in  the  Fermi  coordinates.  The  speed  of  light  is  isotropic  (the  same 
in  all  directions),  but  is  not  homogeneous.  Specifically,  the  speed  of  light  depends  on  Fermi  coordinate  position  and 
time,  but  not  on  direction.  Equation  (83)  gives  the  change  in  the  speed  of  light  from  c,  its  vacuum  value  in  an  inertial 
frame. 

I  have  computed  the  magnitudes  of  these  effects  for  numbers  relevant  to  GPS  satellite  orbits  and  orbital  frequency. 
The  effects  are  all  small  for  GPS,  except  for  the  12  ns  periodicity  for  an  orbiting  clock,  given  in  Eq.  (67),  that  occurs 
at  the  difference  frequency,  —  a;,  between  satellite  orbital  frequency  Ct  and  Earth  orbital  frequency  lo.  This  effect  is 
reminiscent  of  the  periodic  effect  in  GPS  described  in  the  introduction.  However,  the  effects  here  are  all  kinematic, 
since  gravitational  fields  are  neglected.  The  inclusion  of  gravitational  fields  is  expected  to  significantly  modify  the 
results.  In  particular,  from  the  work  of  Ashby  and  Bertotti,  this  12  ns  periodic  effect  is  expected  to  cancel  when 
gravitational  fields  are  included  [2,28].  Work  is  in  progress  to  extend  these  results  to  include  the  effects  of  gravitational 
fields. 

While  the  effects  discussed  in  this  work  are  small  for  the  present  GPS,  these  effects  are  quite  significant  for  future 
navigation  and  time  transfer  systems  that  are  expected  to  have  sub-nanosecond  time  transfer  accuracy  [6-8]. 
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APPENDIX  A:  CONVENTIONS  AND  NOTATION 


I  use  the  notation  and  conventions  of  Synge  [12],  except  that  my  space-time  indices  on  coordinates,  take  values 
i  =  0, 1,2,3.  Specifically,  I  use  the  convention  that  Roman  indices  take  the  values  i  =  0, 1,2,3,  and  Greek  indices 
Q  =  1, 2,3.  I  use  the  summation  convention  for  repeated  indices:  when  the  same  index  appears  in  a  lower  and  upper 
position,  summation  is  implied  over  the  range  of  the  index. 

If  and  x*  H-  dx^  are  two  events  along  the  world  line  of  an  ideal  clock,  then  the  square  of  the  proper  time  interval 
between  these  events  is  dr  =  ds/c,  where  ds  is  given  in  terms  of  the  space-time  metric  as 

ds^  =  —Qij  dx*  dx^  (Al) 

and  I  choose  gij  to  have  the  signature  -1-2.  When  gij  is  diagonalized  at  any  given  space-time  point,  the  elements  can 
take  the  form  of  the  fiat  space  Minkowski  metric  given  by  poo  =  —1,  9cx(3  =  and  poa  =  0. 

In  this  work,  I  take  the  global  inertial  frame  coordinates  to  have  the  Minkowski  metric.  The  observer’s  frame  is 
described  by  Fermi  coordinates,  a  =  0, 1,2,3,  with  line  element 

ds^  =  -G^ij)dX^^Ux^'>  (A2) 
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FIG.  1.  The  observer’s  circular  3-d  path  is  shown  against  the  background  of  the  global  inertial  Minkowski  coordinates  x\ 
The  position  of  the  observer  (origin  of  Fermi  coordinates)  and  the  orientation  of  the  Fermi  coordinate  axes  are  shown 
schematically  at  a  given  time. 
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h  the  initial  tetrad  basis  vectors  Aj,  at  s  =  0  at  point  Pq.  The  Fermi 
s  are  shown  at  point  F. 
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FIG.  3.  The  global  inertial  coordinate  axes  x*  and  Fermi  coordinates  are  shown  in  relation  to  the  coordinate  position 
X,  unit  vector  n,  and  the  plane  X  •  n  =  0,  at  a  given  coordinate  time 


TABLE  L  Numerical  values  of  constants. 


constant 

definition 

value 

a 

Earth  orbital  semi-major  axis 

1.496x10“  m 

U) 

Earth  orbital  angular  velocity 

1.99238x10-^  sec"^ 

b 

GPS  satellite  semi-major  axis 

2.656177x10^  m 

n 

GPS  satellite  angular  velocity 

1.45842x10-^  sec-i 

c 

vacuum  speed  of  light 

2.997924x10®  m/s 

